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Abstract 

Measurements of on-shell ZZ production are described, using data collected 
by DELPHI in 1997 and 1998, at centre-of-mass energies ^Js = 182.6 GeV 
and 188.6 GeV respectively. Results obtained in each of the final states qqqq, 
H~^fi~qq, e'^e'qq, wqq, and ui'l'^l' are presented. The measured 

cross-sections for on-shell ZZ production via the tree-level doubly-resonant 
graphs (NC02) are: 

aNC02 (182.6 GeV) = 0.38 ± 0.18 (stat) ± 0.04 (syst) pb, 

crNC02 (188.6 GeV) = 0.60 ± 0.13 (stat) ± 0.07 (syst) pb. 

They are consistent with the Standard Model expectations of 0.25 pb and 0.65 
pb at each energy. 
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1 Introduction 

The study of doubly resonant production of Z bosons is a relatively new topic. The 
first evidence was accumulated during 1997, when LEP was operated at a centre-of- 
mass energy of 182.6 GeV, corresponding to the threshold for this channel. In this letter 
we present measurements of the production cross-section both from that run and from 
the 1998 run at 188.6 GeV. 

There are several motivations for studying this channel. Firstly, it is necessary to 
check that the observed production rate of any expected physical process seen for the first 
time in a new energy domain is well accounted for by the Standard Model. Deviations 
from predictions could be interpreted as a signal for new physics beyond the Standard 
Model, manifesting itself through anomalous production p|J^, for instance by means of 
anomalous neutral-current triple gauge boson couplings 0. Secondly, ZZ production 
forms an irreducible background to the Higgs search at LEP when the mass of the Higgs 
boson is close to that of the 2' 0. In this context the results obtained also give some 
indication of the reliability of the techniques used in the Higgs search. 

In what follows, the data sets and simulations used are described and the signal def- 
inition which was adopted is discussed. The event selections developed in the six sub- 
channels which were analysed are then presented. Results are given in the form of a 
comparison of the numbers of found and predicted selected events, together with an eval- 
uation of the main systematic effects. Finally, the combination of the sub-channel results 
into overall ZZ cross-sections is described and the overall measurements are compared 
to the Standard Model expectation. 



2 Data samples 

DELPHI took data at centre-of-mass energies of 182.6 GeV in 1997 and 188.6 GeV 
in 1998, with integrated luminosities of 54 pb~^ and 158 pb^^ respectively. A detailed 
description of the detector and a review of its performance can be found in f^J^. The 
detector was not changed in recent years, except for upgrades of the vertex detector 
and the addition of a set of scintillator counters to veto photons in blind regions of the 
electromagnetic calorimetry, at polar angles near 40° and 90° . 

Simulated events were produced with the DELPHI simulation program DELSIM[0 and 
were then passed through the same reconstruction and analysis chain as the data. The 
generation of processes leading to four-fermion final states was done with EXCALIBUR[^], 
relying on JETSET 7.4 for quark fragmentation. GRC4F|11[ was used as a complemen- 



tary generator for four-fermion final states resulting from WeUe processes when cos^e 
> 0.9999. Two fermion processes e^e^ //(+^7) were generated using PYTHIA |T0[, 
e'^e" fi^ fi^ {+n'~f) and e"'"e~ —>■ T~^T^{+n'y) with KDRALZ||12||, and e"'"e~ — * e^e^l+wf) 
with BHWIDE[0|. Two-photon interactions were generated using TWOGAM [1^] and BDK 



3 Signal definition 

The region of phase-space at high di-fermion masses must be isolated to measure the 
ZZ production cross- sect ion. In order to interpret the measurement in terms of the 
tree-level doubly-resonant graphs shown in Figure ^ (referred to as the NC02 graphs) 
the presence of other four-fermion processes in this region must be taken into account. 
For this purpose the signal was defined in the simulation by requiring that the generated 



2 



masses of the two appropriate pairings of final-state fermions be within 10 GeV/c^ of the 
nominal Z mass. This choice of mass window maximized the sensitivity to the NC02 
graphs while minimizing contributions from other four-fermion processes. Events with 
the correct flavour composition but which fell outside this generator-level mass window 
were considered as background. A scaling factor R, 

^total 

j^^^NO)^ (1) 
-.window ' ' 

was then calculated at generator level to enable conversion of the measured total four- 
fermion cross-section within the mass window, cr^™^°™, into the total NC02 cross-section, 
(j^co2- The scaling factors obtained for each channel are shown in table |l[ Most of the 
values are close to unity, confirming that the defined region is dominated by on-shell 
ZZ production. Several sources of bias which could result from this procedure were 
investigated. The most relevant among these were: 

• In the channel e^e~qq, the scaling factors are smaller than unity. This arises be- 
cause of significant contributions from the single- resonant process e^e~ Z. Since the 
electrons from this process tend to be peaked in the forward directions and since the 
efficiencies to identify electrons are reduced in these regions, there is a bias for this 
particular channel from taking a scaling factor averaged over the full solid angle. 
The magnitude of this bias was estimated by computing the expected cross-sections 
using this procedure in the barrel and forward regions separately. A 2.75% correction 
was derived and incorporated into the scaling factor quoted in table 

• In the channel l^l'l^l^ the factor is significantly smaller than unity. This arises 
because of the intrinsic ambiguity in pairing, existing in cases with four leptons of 
the same family, for which there are significant contributions in the signal window 
from Z7* processes. As an additional cross-check for this channel the measurement 
was also repeated using an additive correction procedure (see the corresponding 
section), yielding fully consistent results. 



Final State 


183 GeV 


189 GeV 


qqqq 


1.26 


1.07 


vvqq 


1.14 


1.05 


fi^fi-qq 


1.19 


1.05 


e~^e~qq 


0.83 


0.88 


l+l-l+l- 


0.46 


0.59 




1.00 


1.00 



Table 1: Scaling factors R computed at generator level to convert the measured total 
four-fermion cross-section within the signal defining mass window \Mjj — Mz\ < 10 
GeV/c^, aJi'"^°'", into the total NC02 cross-section, a'^^^''^^- 



4 Four jets 

The ZZ qqqq process represents 49% of the ZZ final states and results typically in 
events with four or more jets. The principal backgrounds arise from WW and qq(7) pro- 
cesses, which can lead to similar multijet topologies. The main ingredients used to isolate 
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the signal were the 6-tagging of the jets, topological information quantifying their sep- 
aration, and reconstructed di-jet masses. A pre-selection was first applied to isolate 



hadronic events with at least four reconstructed jets, as in Four and five jet events 
were treated separately throughout the analysis^. 

In order to optimally use the information from the 6-tagging, the ratio of the distri- 
butions of the combined 6-tagging variable [O] for b and non-6 jets were parametrized in 
three different angular regions, using simulated Z decays. Using these parametrisations 
and the measured value of the combined 6-tagging variable for each jet, the probability 
that it originated from the fragmentation of a 6 or a non-6 quark was evaluated. 

In order to fully exploit the information from the mass reconstruction, all possible 
di-jet pairings were used, following the method outlined in |^6|. For each of the physics 
hypotheses ZZ , WW and qq(7) , the expected probability distributions for the different 
mass combinations in an event were expressed in analytical form, by means of products 
of Breit-Wigners and phase-space factors for the correctly paired mass combinations in 
WW and ZZ events, and of flat spectra for all the other cases. Using the event-by- 
event jet errors, measured probability distributions were obtained, as functions of 



the possible di-jet masses, by kinematic fits |T8[ requiring four-momentum conservation 
and equality of each di-jet mass with the tested value. These two-dimensional measure- 
ment probabilities were then convoluted with the expected distributions to quantify the 
compatibihty of each pairing with the ZZ WW and qq(7) hypotheses from the mass 
information alone. 

Finally the topological variable Emin ■ amin was defined, where Emin is the minimum 
jet-energy and amin the minimum opening angle between any jet pair. The ratio of two- 
fermion and four-fermion events was parametrized as a function of this variable using 
simulated data for four- and five-jet events separately. 

Looking for all possible hadronic ZZ final states and using the predicted SM cross- 
sections and branching ratios into the different quark configurations, the 6-like probability 
per jet, the topological information per event and the mass information per pairing, a 
combined variable quantifying the compatibility of any event with the ZZ hypothesis 
was constructed. This combined variable was shown to behave as a genuine probability 
for the ZZ hypothesis, and provided high discriminating power. The distribution of this 
variable at 188.6 GeV is shown in Figure ^ Results obtained by cutting on this variable 
are presented in table ^j. 



182.6 GeV 


data 


Total MC 


ZZ 


ZY 


WW 


qq(7) 


After pre-selection 
At maximum e-p 


535 
2 


510.8 
1.68 


6.73 
1.01 


7.06 
0.10 


336.2 
0.11 


160.8 
0.46 




188.6 GeV 


data 


Total MC 


ZZ 


Z-f* 


WW 


qq(7) 


After pre-selection 
At maximum e-p 
At S/B = 2.92 


1538 
42 
7 


1536.0 
43.06 
7.69 


41.8 
15.54 
5.73 


17.7 
1.41 
0.38 


1039.6 
16.01 
0.28 


437.0 
10.10 
1.30 



Table 2: The observed and expected number of selected events after the pre-selection 
at 182.6 GeV and at 188.6 GeV, for integrated luminosities of 53 pb^^ and 158 pb^^ 
respectively. The two last lines give the number of events left at the maximum of efficiency 
X purity (e-p) and at a signal-to-background ratio, S/B, of about 3. 



Events with more than five reconstructed jets were forced into a five-jet configuration. 
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Systematic effects resulting from uncertainties in signal efficiency, in the conversion 
factor to translate the results to a NC02 cross-section and in the predicted background 
level were studied. The largest errors came from uncertainties in the modeling of multijet 
qq(7) processes with b quarks composing the main remaining background, and from the 
sensitivity of the 6-tagging procedure to the jet multiplicity. The dominant qq(7) mod- 
eling uncertainty came from the limited precision with which the g ^ bb rate is known, 
and to a lesser extent from uncertainties in the amount of reduction in gluon radiation 
off b quarks (also known as the dead-cone effect). These uncertainties, as well as biases 
in the treatments provided by the generator, have been studied |]19[ by comparing results 
of analytic calculations with dedicated measurements at LEP-1 and with predictions of 
the generators used in the simulation. Uncertainties in selection efficiencies related to the 
6-tagging procedure were studied by comparing efficiencies for two and four jet events at 
LEP-1 energies. By propagating effects from these sources into the analysis, a combined 
systematic uncertainty of around 6% was estimated relative to the expected cross-section. 

5 Jets and a pair of isolated leptons 

The decay modes e~^e^ f^^fJ' QQ and e+e~ —>■ e^e~qq represent 9% of the ZZ final 
states. Events with r+r^ pairs were not considered. The two final state leptons are 
typically well isolated from all other particles. This can be used to select such events 
with high efficiency in both the muon and electron channels. Events were selected initially 
without explicit cuts on the masses of the final state fermion pairs in order to analyse 
simultaneously ZZ ^ Z7* events and contributions from other possible diagrams leading 
to final states such as Ze^e~ or t-channel 7* exchange with Z/7-strahlung. Mass cuts 
were then applied to isolate the ZZ contribution. 

A loose hadronic pre-selection was first applied, requiring that events have at least 
7 charged particles and a charged energy above 0.30 a/s. To suppress the radiative 
return to the Z boson, events were rejected if a photon with energy more than 60 GeV 
was found. The selection procedures then proceeded in a closely similar way for both 
H~^IJ'~qq and e^e~qq channels. In order to maximize the lepton identification efficiency, 
any charged particle with a momentum exceeding 5 GeV/c was considered as a possible 
lepton candidate around which nearby photons, if present, could be clustered. This 
was found to be necessary to improve the energy evaluation in the presence of final 
state radiation, and, in the case of electrons, bremsstrahlung. In the case of the e^e~qq 
channel, photons with energy between 20 GeV and 60 GeV were also considered as 
electron candidates, to recover events in which the electron track was not reconstructed. 

Events with at least two lepton candidates of the same flavour, opposite charge and 
invariant mass exceeding 2 GeV/c^ were then selected. All particles except the lepton 
candidates were clustered into jets and a kinematic fit requiring four-momentum conser- 
vation was applied, correcting appropriately the errors on lepton energies in cases where 
photons had been added by the clustering procedure. 

At least one of the two lepton candidates was required to satisfy strong lepton iden- 
tification criteria, while softer requirements were specified for the second. Muons were 
considered as strongly identified if selected by the standard DELPHI muon identification 
package |^, based mainly on finding associated hits in the muon chambers. For soft 
muon identification only a set of kinematic and calorimetric criteria were used. Electrons 
were considered as strongly identified when the energy deposited in the electromagnetic 
calorimeter exceeded 60% of the cluster energy or 15 GeV and when the energy deposited 
in the hadron calorimeter was reasonably limited. For soft electron identification only 
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requirements on the momentum of the charged particle in the cluster and on the energy 
deposited in the hadron calorimeter were used. Moreover electron candidates originat- 
ing from applying the clustering procedure around a photon were considered as softly 
identified. 

Two discriminating variables were then defined for final event selection: p™"^ the 
lesser of the transverse momenta of the lepton candidates with respect to their nearest 
jet and the per degree of freedom of the kinematic fit. The distribution of the mass 
of one fermion pair (/+/^ or qq) when the mass of the second pair is within 15 GeV/c^ of 
the nominal Z mass is shown in figure ^,b. The distribution of the sum of masses of two 
fermion pairs is shown in figure The observed distributions are in reasonable agree- 
ment with the predictions from simulation. To select on-shell ZZ production, cuts were 
placed simultaneously on the masses of the pair, on the remaining hadron system, 
and on their sum, taking into account in the performance optimization the different mass 
resolutions of these final states and the presence of the single-resonant process, e^e~ Z in 
the case of the e^e~qq channel. The observed and predicted numbers of selected events 
are shown in table |^. The background is divided into two parts: the first one comes from 
l^l~qq events outside the generation- level signal window, and the second one comes from 
other processes, principally W~^W~, other ZZ decays and (in the case of e^e~qq) qq{+'^) 
production. For the fi^fi~qq channel, the efficiencies of the selection were 0.89 ± 0.02 
and 0.86 ± 0.01 at 182.6 and 188.6 GeV respectively. For the e^e~qq channel, they were 
0.73 ± 0.03 and 0.72 ± 0.02 respectively. 



E 


channel 


data 


Total MC 


MC signal 


MC l+l-qq 


MC other 


GeV 










backgr. 


backgr. 


182.6 


fi^fi'qq 


3 


0.52 ±0.03 


0.48 ±0.03 


0.04 ±0.01 


0.00 ±0.01 




e~^e~qq 





0.69 ±0.07 


0.55 ±0.05 


0.10 ±0.02 


0.04 ±0.04 




Total 


3 


1.21 ±0.08 


1.03 ±0.06 


0.14 ±0.02 


0.04 ±0.04 


188.6 


^l^H~qq 


5 


4.15 ±0.11 


3.92 ±0.10 


0.19 ±0.02 


0.04 ±0.03 




e~^e~qq 


3 


4.04 ±0.15 


3.55 ±0.13 


0.35 ±0.04 


0.14 ±0.06 




Total 


8 


8.19 ±0.19 


7.47 ±0.16 


0.54 ±0.04 


0.18 ±0.07 


All 


fi+fi-qq 


8 


4.67 ±0.11 


4.40 ±0.11 


0.23 ±0.02 


0.04 ±0.03 




e~^e~qq 


3 


4.73 ±0.17 


4.10 ±0.14 


0.45 ±0.05 


0.18 ±0.07 




Total 


11 


9.40 ±0.20 


8.50 ±0.18 


0.68 ±0.05 


0.22 ±0.08 



Table 3: The observed and expected number of selected l^l^qq events - candidates to 
the on-shell ZZ production. The errors shown are due to the simulation statistics only. 
The expected numbers of the signal events are for the SM cross-section of the on-shell 
ZZ production. 

Several sources of systematic errors were investigated. Uncertainties in the lepton 
identification were estimated comparing semileptonic W^W~ events selected in data and 
simulation using the strong lepton identification criteria. Uncertainties in signal efficien- 
cies from the description of the kinematic observables used were evaluated comparing 
the Pt and distributions in data and simulation for all l^l~qq events selected without 
mass cuts. Corresponding uncertainties in background levels were evaluated by compar- 
ing samples of events selected in data and in simulation requiring both isolated tracks 
not to be identified as leptons, while maintaining all other criteria. Finally, uncertainties 
in the background level in the e^e~qq channel from fake electrons were studied with Z7 
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events selected in data and in simulation with purely kinematic criteria. Propagating 
the found differences to the final stage of the analysis yielded a combined systematic 
error on the efficiency to select l^l^'qq events of ±3.0%. The uncertainty resulting for the 
background level amounted to about ±15%. 

6 Jets and missing energy 

The decay mode vvqq represents 28% of the ZZ final states. Its signature is a pair of 
jets relatively acoplanar with the beam and with visible and recoil masses compatible with 
the Z mass. The most difficult backgrounds arise from single resonant Wev^ processes, 
from WW processes where one of the W decays into rvr-, and from events, accompanied 
or not by isolated photons escaping detection, in which one or two of the jets were badly 
reconstructed. 

A pre-selection was first applied to remove the bulk of the background. Multihadronic 
annihilation events were selected by requiring that the number of charged particles be 
larger than 8, that the track of at least one charged particle, with a transverse momentum 
larger than 1.5 GeV/c, extrapolate back to within 200 /xm of the primary vertex in the 
plane transverse to the beam axis, that the total charged energy of the event exceed 10% 
of the centre-of-mass energy and that its raw visible mass be in the range 81 ± 22 GeV/c^. 
To reject radiative returns to the Z with energetic photons emitted in blind regions of the 
electromagnetic calorimetry (at polar angles near 40° and 90° ) signals from dedicated 
scintillator counters were used. 

A combined discriminant variable was then constructed using the Iterative Discrim- 
inant Analysis program (IDA) |2y] to calculate a second order polynomial from seven 
event variables, selected based on their discriminating power and independence: 

• The minimum polar angle defining a cone in the positive and negative beam direc- 
tions containing 15% of the total visible energy. 

• The logarithm of the acoplanarity scaled by the sine of the minimum polar angle 
between a jet direction and the beam axis. 

• The total reconstructed energy. 

• The maximum transverse momentum between any particle and a jet. 

• The energy of the more energetic among the reconstructed jets. 

• The thrust, computed in the rest frame of the visible system. 

• The acoUinearity of the jets. 

Moreover, to concentrate on the signal region when optimizing the second-order dis- 
criminant function, very loose cuts were applied on these variables, at values correspond- 
ing to the tails of the signal distributions and removing each about one percent of it. 
Finally, the total energy of hypothetical photons escaping in the beam direction was es- 
timated from the two jet directions and was required to be less than 85% of the value 
expected for radiative returns to the Z. 

The comparison of selected data and simulation rates for the signal and background 
components is shown in Figure ^ (for the 188.6 GeV data) and in table ^ The ZZ 
cross-section at each energy was obtained from a binned maximum likelihood fit of this 
discriminant output with the ZZ signal contribution as the only free parameter. 

At both energies more events than predicted were observed before the final selections. 
This resulted mainly from uncertainties in the description of the energy reconstruction 
of the gg(7) background. Unlike the other ZZ decay channels studied, in which the 
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loo yjcG V 


uaxa 


loxai ivn^ 


Signal 


4-1 




qql^l 


77 


Xi.1 LCI IJl C oClCl^ IjlWli 




2286+1 n 


3 5+n 1 

(J . O -i—\J . ± 


1 34+1 


2005+4 


6 9+0 2 

\J . U -i—\J . ^ 


1 37+9 


After tail cuts 


235 


209±2 


3.0±0.1 


37+0.6 


163+1 


1.0+0.1 


5.1+1.6 


After final DA cut 


22 


21.9±0.45 


2.4±0.1 


8.8+0.3 


10.6+0.4 


0.07+0.03 





189 GeV 


Data 


Total MC 


Signal 


CC 4-f 






77 


After pre-selection 


6723 


6206±20 


21.5±0.3 


344+2 


5344+7 


17.1+0.5 


479+18 


After tail cuts 


824 


710±6 


20.4±0.3 


158+1 


500+3 


3.3+0.1 


28.5+4.6 


After final DA cut 


112 


106±2 


17.0±0.3 


47.5+1.1 


39.2+0.7 


0.5+0.05 


1.5+1.0 



Table 4: vvqq channel: Data and simulation rates after different steps of the analysis 
(see the description in the text). The uncertainties result from the limited simulation 
statistics. The CC 4-f background refers to charged-current four-fermion processes such 
as W+W-. 



full event energy is reconstructed most of the time, the selection of ZZ — ^ uuqq events 
exploits the large missing energy characteristic of this channel, and is hence sensitive to 
the description of the low energy tail in the reconstruction of background processes such 
as qq{'y)- To study such effects and evaluate their propagation to the final steps of the 
analysis, large statistics samples of Z events, collected in the same conditions as the high 
energy data, were compared to the simulation to estimate corrections to the particle flow, 
in bins of momentum, polar angle and particle type. These consisted mainly of changes 
in multiplicities, to account for observed efficiency losses and possible duplication effects 
in the pattern recognition. The corrections, typically of one to a few percent in the 
barrel and endcaps respectively, were largest in the case of neutral particles and for high 
reconstructed momenta, and resulted in significantly improved agreement between data 
and simulation for energy flow observables such as the total charged and neutral energies, 
or the visible mass, including in the tails of these distributions. This procedure was then 
applied to all simulated high energy samples. The distribution of the discriminating 
variable obtained at 188.6 GeV after applying these corrections is shown in Figure ^3. 
The main effect was to increase the gg(7) background, by up to 15%. The ZZ cross- 
section fit was then repeated using this modified version of the simulation and differences 
were used as a conservative measure of the impact of uncertainties in the description 
of the low energy tail of the reconstruction. To reduce effects, regions where the gg(7) 
background dominated were removed by a final cut on the discriminating variable at zero 
(indicated by an arrow in Figure ^). Above this cut, 22 events remained at 182.6 GeV 
and 112 events at 188.6 GeV. The systematic shift in the cross-section obtained amounted 
to 30% relative to the expected value, and was used to represent the systematic error at 
both energies. 

7 Four-lepton final states 

About 1% of all ZZ events lead to the I'^l^l^l^ final state. The event topology is 
clean and the only significant background comes from non-resonant e^e~l^l~ production. 
Events were selected if they contained between 4 and 8 charged particles, accompanied 
by at most 10 neutral particles, irrespective of particle identification. In order to take 
into account final state radiation and bremsstrahlung effects for candidate electrons, 
the momenta of the charged particles were rescaled if the measured sum of energies of 
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electromagnetic clusters in a narrow cone around the track direction was larger than the 
energy inferred from the track momentum measurement. The total invariant mass of the 
charged particles had to be greater than 50 GeV/c^, and the minimum invariant mass 
after discarding any one of the charged particles larger than 20 GeV/c^. All combinations 
of four charged particles with total charge zero were then examined, and a combination 
was selected if: 

• all four tracks had their impact parameters at the interaction point smaller than 3.0 
and 0.5 cm, respectively in the projections containing the beam axis and perpendic- 
ular to it, and polar angles between 10° and 170°, 

• at least three of the four charged particles had momenta greater than 5 GeV/c, and 
the least energetic particle a momentum greater than 2 GeV/c, 

• a system of two oppositely charged particles was found with both their invariant and 
recoil masses within 10 GeV/c^ of the Z boson mass and having the same flavour 
(when both were identified as either electrons or muons), 

• the two particles complementary to this system were separated by at least 90° from 
each other, 

• the invariant mass of all pairs of oppositely charged particles in the event exceeded 
2 GeV/c2. 

At 182.6 GeV, no event was found, where predicted signal and background of respec- 
tively 0.098 ± 0.005 and 0.041 ± 0.103 events are expected. At 188.6 GeV one event 
was observed, where predicted signal and background of respectively 0.683 ± 0.022 and 
0.076 ± 0.043 events are expected. The efficiencies of the selection were respectively 
0.31 ± 0.03 and 0.41 ± 0.02 at each energy. 

To study potential systematic effects arising from the method based on mass windows 
used to disentangle the NC02 contribution from other non-resonant four fermion pro- 
cesses, the results were checked with a different signal definition, consisting of generating 
two samples, the first using NC02 graphs only and the second using all other 

tree-level graphs, with NC02 excluded. The first sample was used to define the signal 
and the second was used to estimate the background from other four-fermion processes. 
Results obtained in this way were fully consistent with those obtained using the mass 
window method. 

8 Two isolated leptons with missing energy 

Events from ZZ — > i/P/j,'^ ,i'i?e~^e~ processes are characterized by two relatively 
acoUinear charged leptons of the same flavour, with both invariant and recoil masses 
close to that of the Z, and by large missing energy. Although it has different production 
kinematics, the WW process also contributes to these final states with a large cross- 
section. A significant fraction of the corresponding events have exactly the same features 
and constitute a dominant background. 

To ensure good reconstruction, tracks were required to have impact parameters at the 
interaction point smaller than 3.0 and 0.5 cm, respectively in the projections containing 
the beam axis and perpendicular to it, and polar angles between 20° and 160°. As for the 
channel, final state radiation and bremsstrahlung effects for candidate electrons 
were taken into account by rescaling the momenta of the charged particles if the measured 
sum of energies of electromagnetic clusters in a narrow cone around the track direction 
was larger than the energy inferred from the track momentum measurement. Events with 
two particles identified as e'^e~ or were selected if their total energy was less than 
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60% of that of the centre-of-mass, if the angle between them was in the range 140°-170°, 
if the polar angle of the missing momentum vector was between 25° and 155°, and if the 
reconstructed invariant masses satisfied: 

min{| Mz - m{l+l-) \ , \ Mz - rrirecoiiil^l') |} < 4 GeV/c^, and 

max{| Mz - m{l+l-) \ , \ Mz - nirecoii{l-^l-) |} < 8 GeV/c^, 
where mrecoii{l^l~) is the invariant mass recoiling against the l^l" pair. 

At 182.6 GeV, no event was found, with predicted signal and background of 0.109 ± 
0.007 and 0.348±0.155 events respectively. At 188.6 GeV, two events were observed, with 
predicted signal and background of 0.839 ± 0.028 and 1.025 ± 0.263 events respectively. 
The efficiencies of the selection were 0.26±0.02 and 0.30±0.02 at each energy respectively. 

The main systematic effects in this channel were from uncertainties in the lepton 
identification efficiencies, both of which are at the few per cent level, and to a lesser 
extent from uncertainties in the WW cross-section used. In view of the limited statistics 
in this channel, these effects were neglected. 

9 Combined NC02 cross-section 

The cross-sections in the qqqq and vvqq channels were determined individually by 
means of binned maximum likelihood fits to the distributions of the combined variable 
defined in each of the corresponding analyses. To derive a combined value for the NC02 
cross-section from all the six channels analysed, a global likelihood was constructed by 
combining the likelihoods from the fits performed in these two channels with the Poisson 
probabilities for observing the number of events seen in the four other channels {fJ'~^fi~qq, 
e~^e~qq, and uul^l'), given the predicted numbers. This global likelihood 

was then maximised with respect to variations of the value of the NC02 cross-section, 
assuming branching ratios of the Z fixed to those expected in the Standard Model, and 
taking into account the scaling factors listed in table |l]. 

The impact of the systematic uncertainties affecting the different final states was 
studied by repeating the fits with modified assumptions on backgrounds and efficiencies 
(see the corresponding sections). Statistical errors exceeded the estimated systematic 
ones in each of the channels taken separately. To estimate the total systematic error 
on the combined result, the different uncertainties were assumed, conservatively, to be 
correlated, and their effects were added. This resulted in a total error of ±11% relative 
to the expected value, dominated by the uncertainty in the uuqq channel. 

The values for the NC02 cross-sections obtained were: 

aNC02 (182.6 GeV) = 0.38 ± 0.18 {stat) ± 0.04 (syst) pb, 

(T7VC02 (188.6 GeV) = 0.60 ± 0.13 (stat) ± 0.07 (syst) pb. 

These values are consistent with the Standard Model expectations of 0.25 pb and 0.65 
pb at each energy, respectively, predicted by EXCALIBUR. The agreement is illustrated in 
Figure ^. Measurements of on-shell ZZ production by the three other LEP collaborations 
can be found in 0,^^. 
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Figure 1: The Feynman graphs for on-shell ZZ production (referred to as the NC02 
graphs) 
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Figure 2: Distribution of the ZZ probability for all events at 188.6 GeV in the four-jet 
channel 
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Figure 3: a) the distribution of the mass of the pair when the mass of the hadron 
system is within 15 GeV/c^ of the nominal Z mass; b) the distribution of the mass of 
the hadron system when the mass of the pair is within 15 GeV/c^ of the nominal 
Z mass; c) the distribution of the sum of the masses of the pair and of the hadron 
system. The points are the data taken at 182.6 GeV and at 188.6 GeV, and the histogram 
is the simulation prediction. The contribution from the signal, as defined in section 3, 
is indicated by the empty histogram. The contributions from backgrounds are indicated 
with light grey when arising from from l^l~qq final states, and with dark filling otherwise. 
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Figure 4: vvqq channel: The second-order discriminant for the data at 188.6 GeV, for 
the expected signal (light) and for the backgrounds from l^l~qq and Z7* (light grey), 
qqij) (dark grey) and Wev^ and WW (dark, CC 4-f in the legend) processes are shown, 
using the nominal (plot a) and corrected (plot b) simulation, following the procedure 
described in the text. The main effect from the correction procedure was to increase the 
99(7) background component. The arrow indicates the cut used to remove the part of the 
distribution where the gg(7) background dominates, in order to reduce the systematic 
error as described in the text. 
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Figure 5: Combined NC02 cross-sections measured from data collected in 1997 and 1998. 
The solid curve was computed using the YFSZZ Monte-Carlo Values obtained with 

the EXCALIBUR Monte-Carlo were within ±1 % of those obtained with YFSZZ. 



